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ABSTRACT 



> 

\^ , We present two-dimensional (2D) mapping of the stellar velocity field within the 

' inner 5" of six nearby active galaxies, using spectra obtained with the Integral Field 

I Unit of the GMOS instrument at the Gemini North telescope. The sampling of the 

, observations is 0'.'2, corresponding at the galaxies to spatial extents ranging from 10 

to 30 pc. The spatial resolution range from 20 to about 180 pc, and the observed 
field of view covers a few hundred parsecs around the nuclei. The Calcium II triplet 
r"| . absorption features at « 8500A were used to measure the stellar radial velocities 

{Dh' and velocity dispersions. The radial velocity fields arc dominated by rotation in all 

galaxies. A simple kinematical model assuming a purely rotating system with circular 
^ • orbits in a plane was fitted to the radial velocity data. The turnover of the rotation 

^ I curve is at only « 50 pc for NGC4051 and between 200 and 700 pc for the other 5 

. galaxies. The velocity dispersion (a) maps show the largest values (100 ^ ct ^ 150 

kms~^) at the centre. In the cases of NGC 2273 and NGC 3227, there is a decrease to 
a « 70 — 80 km s~^ at « 200 — 300 pc from the nucleus, delineating partial rings of low 
a values. A similar broken ring seems to be present at ~ 400 pc from the nucleus also 
' in NGC 4593. We interpret these low a rings as traces of recently formed stars that 

5^ , partially keep the cold kinematics of the original gas from which they have formed. In 

NGC 3516 there is a decrease of cr outwards with the steepest gradient observed along 
the direction of the galaxy major axis, where a reaches « 80 — 90 kms^^ at « 400 pc 
from the nucleus. 

The main novelty of the present work is the unprecedented spatial resolution 
reached by a 2D study of stellar kinematics of Seyfert galaxies using an IFU. The few 
similar IFU studies available in the literature for Seyfert galaxies have a much poorer 
spatial resolution and/or are restricted to the study of emission line kinematics. 

Key words: galaxies: active - galaxies: Seyfert - galaxies: nuclei - galaxies: starburst 
- galaxies: kinematics and dynamics - stellar dynamics 



1 INTRODUCTION 



It is now widely accepted that active galactic nu- 
clei (AGN) are powered by the accretion of material onto 
a central supermassive black-hole (SMBH). The present 
paradigm for the evolution of galaxies is that all galax- 
ies which form bulges also form a SMBH at the ir centres 
JPerrarese fc Merritd lioool : iGebhardt et al.ll2000l) . In this 



scenario, the active galaxies are those in which the SMBH 
is presently accreting material from its surroundings. 

A problem still under investigation is the mechanisms 
by which the material is dragged to the nuclear regions 
to feed the AGN, as the galactic gas must loose almost 
all of its angular momentum in order to reach the cen- 
tral few parsecs. Signatures that such feeding is occurring 
include the frequent occurrence of recent star formation 
around AGN, whi ch implies the existence of a gas reservoir 
close to the AGN jSchmitt et alJll999l : ICid Fernandes et alJ 
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| 200ll: IStorchi-Bergmann et al.l I2OO0I: iBoisson et alJ I2OO0I: 
ICid Fernandes et alJ l2005l : IStorchi-Bergmann et al.l l20m 

Besides the signatures of young stars observed in spectra of 
AGNjkiloparsec scale kinematic studies (Nelson fc Whittle 
Il995l . ll99a) have suggested that Seyfert galaxies have lower 
mean mass-to-light ratios than normal galaxies, which also 
could be due to a youn ger near-nuclear stellar population 
JOlivaet al.lll995l. 'l999\ 

The presence of young stars near the nucleus seems also 
to be the explanation for the results of recent studies in- 
vestigating the stellar kinematics on sca les of hundred of 
parsecs |Ems ellem et al. 2001; Garci'a-L orenzo et al.llT999l : 
iMarauez et alJl2003ll . These studies have found a central 
drop in the stellar velocity dispe rsion (hereafter g-dro p) in 
a few galaxies with Seyfert nuclei. IWozniak et alJ (|20o3) pre- 
sented simulations showing that a new generation of stars 
formed at the centre of the galaxy from cold material would 
create such a drop which would remain visible for hun- 
dreds of Myrs. Nevertheless, a decrease in a towards the 
nucleus has be en recently o b served also in late-type non- 
active galaxies iGanda et al.l i200(i) . supporting its link to 
recent star-formation, but not necessarily to nuclear activ- 
ity in galaxies. 

Most studies available in the literature on kinematics of 
Seyfert galaxies are based on long-slit observations, which 
are restricted to only one axis of the host galaxies. In order 
to properly probe the galactic gravitational potential, as well 
as to investigate the nature and extent of the cr-drops, it 
is necessary to cover a 2D region, what we have done in 
the present work using the Integral Field Unit (IFU) of the 
Gemini Multi-object Spectrograph (GMOS). 

The power of IFU observations to map the large scale 
kinemat ics of galaxies has been evidenced in recent stud- 
ies fe.g. iGanda et al."200(fl. For AGN, only few such stud- 
ies are available (e.g. EmscUcm eTaDbOOg for NGC1068). 
Most AGN kinematic studies have focused on gas emission- 
lines, as it is much easier to obtain kinematical data from 
emission lines than from stellar absorption features. Never- 
theless, the gas in the nuclear region of AGNs is subject to 
non-gravitational effects, such as winds and jets, and in or- 
der to probe the gravitational potential, it is necessary to 
measure the stellar kinematics. 

In this work, we use the IFU-GMOS to map the stel- 
lar kinematics in the inner few hundred parsecs of 6 of the 
closest Seyfert galaxies. The combination of the spatial res- 
olution at the galaxies reached by the observations (down to 
a few tens of parsecs), 2D coverage, stellar kinematics, and 
a Seyfert sample is unique in the literature. 

This paper is organized as follows. In Section 2 we 
present the criteria used to select the sample and discuss 
the relevant information for each galaxy. In Section 3 we 
describe the observations and reductions, in section 4 we 
present the methods used to analyze the data and in Sec- 
tion 5 we report and discuss our results. In section 6 we 
present a summary of the results and our conclusions. 



2 SAMPLE 

We chose to map the stellar kinematics using the Ca 
II triplet absorption feature around 8500A (hereafter Ca-T) 



due to the fact that this spectral region is not much affected 
by emission lines and the continuum of the active nucleus, 
and thus clearer measurements of the stellar kinematics can 
be made. 

The sample galaxies were selected as the closest Seyfert 
galaxies for which previou s spectra in the Ca-T reg ion were 
available in the literature jNelson fc WhittlelllQQ^ . so that 
we could check, in advance, if the absorption features to 
be used in the stellar kinematic measurements were clearly 
present in the spectra. In this work, we present the results 
for six Seyfert galaxies, whose properties are listed in Table 
□ 

Besides selecting the galaxies on the basis of their prox- 
imity and detectability of the stellar Ca-T, we also looked for 
the presence of strong [S ill] A9069 line emission, as this line 
could be included in the same observational set up used to 
observe the Ca-T. The results based on the measurements 
of this emission-line, such as intensity maps and gaseous 
kinematics, will be presented in a forthcoming paper. 

We have also observed kinematic standard stars, to be 
used as references in the measurement of the radial velocities 
and velocity dispersions of the sample galaxies. 

2.1 NGC2273 

F rom a 6 cm radio map, lUlvestad fc Wilson! 
il984l) found the nucleus to be resolved into two 
distinct components sepa, r a,ted by 0'.'9 along E-W. 
iFerruit. Wilson fc MulchaevI l)200nh . using HST, resolved 
the central ovoid of the galaxy, which shows two arc-like 
structures forming a partial nuclear ring, with semi-major 
axis of ~ 2'.'5, observed both in a broad-band and in a 
[Nii]-|-Hq emission- line image. They suggest that emission 
in the ring, also weakly observed in [O ill] , comes from H 11 
regions. 

The partial ring is also detected in a color map 
log (F547M/F791W) which presents lower values to NW 
suggesting this is the near side of the ring. The analysis 
of the [Qui] and [Nll]-|-Ha emission line images reveals a 
jet-like structure which extends by 2" to E of the nucleus 
and is aligned with the radio structures observed in 6 cm. 

Subsequent investigation (|Erwin fc SDarkjl200,^ based 
on HST images and colour maps led the authors to classify 
the nuclear ring as a star-forming ring. They also classified 
the structure internal to the ring as a luminous two-armed 
blue spiral in opposition to a previous bar-like interpretation 
jMulchaev fc ReganllT997h . 

The large H — K colour index and high 10 /im lu- 
minosity, extended by more than 5" also supports the 
presence of a circumnu clear starburst foevereux .19891 : 
lYoung fc DevereuJl99 j) . 

2.2 NGC3227 

This galaxy forms an interacting pair with th e ellip- 
tical galaxy NGC 3226. IConzalez Delgado fc PereS ([1993) 
reported traces of young stellar popu lation in the nuclear 
optical spectrum. iMundell et alJ l)l995r) studied the H l emis- 
sion and found that the galaxy disk is inclined by 56°with 
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major axis at PA= 158° and contains a n Hi mass of 5 . 7 x 10* ^ 
Mq. H2 maps have been obtained bv lQuillen et alJ (Il999l) 
from HST/NICMOS images and show elongated emission 
along PA~ 100°.|M eixn cr et al. (1990} obtained ^^CO (1-0) 
maps showing a nuclear double peak aligned roughly E-W 
separated by 2" and extended emission running from SE to 
NW at 30° from the major axis of the galaxy (the direction 
pointing t o the companion galaxy NGC 3226) . To explain 
these maps lf^rnandez et al.i (1999) proposed a gas "disk" of 
100 pc diameter with a major axis along PA~ 100°. Better 
res olution "CO (1-0) and ^^CO (2-1) m aps were obtained 
by ISchinnerer. Eckart fc TacconJ i2000l) who find that the 
central region is resolved into an uneven ring-like structure 
with radius ~ 3" and the E part six times brighter than the 
W part. These maps also show the SE-NW extended struc- 
tures and a nuclear bar that connects the ring to the external 
NW component. The kinematical analysis led them to con- 
clude that a warped gas disk provides a better description 
of the observed gas motions than a bar. 



2.3 NGC 3516 

This was the first Seyfert galaxy with detected 
line variability jAndrillat fc So uffrin|[l96al . More recently 
IWeaver Gelbord^^^aoobril20^i) obs erved variability in 
both the X-ray continuum and Ka line. 'Rega n fc MulchaevI 
11999..) . using HST images showed that a single dust spiral 
patter n dominates the nuc lear morphology. Using J-band 
images |Quillgi^^L|_ IM^lL proposed that this galaxy has 



Table 1. Sample galaxies (column 1), morphological/Seyfert type 
(columns 2/3), adopted spatial scale (column 4) and seeing mea- 



two bars. iNaear et al 



1119991) found a radio jet along P.A. 
~ 10° which correspo nds to the near side o f the galaxy, 
according to Malka n. Goriian fc TamI lll998il . so this jet 
would be projected against the near side of the galaxy. 
Hq -I- [Nil] and [Olll] images show a "Z" -shaped cir- 
cumnuclear emission extended by ~ 20" from SW to NE 
JPoggelll989l: iMivaii. Wilson fc Perez-Fournonlll992h . The 
gas velocities inside this structure cannot be reproduced 
by rotation, but can be explained by a nuclear outflow 
jMulchaev et alJ Il992l: IVeiUeux. TuUv. fc Bland-HawthornI 
l^9^) ~Th e galactic disc does not contai n Hll regions, ac- 
cor dinp^ to [Go nzalez D ejgado et all 1^93). 

lArribas et al.. (.1997ll presented 2D stellar kinematics 
over a similar field to ours but with poorer spatial reso- 
lution. Details about their data and a comparison to ours 
are discussed in Sec. |K1 



2.4 NGC 4051 

Using [O III] emission-line archival HST images, 
ISchmitt fc Kinney (1996) found an unresolved nuclear 
source and a low surface brightness component extending 
l'.'2 from the nucleus along PA = 100°, approximately the 
same orientation as that connecting the two radio compo- 
nents at 6 cm detected bv lUlvestad fc WilsonI ||1984D . which 
are separ ated by 0'.'4. In radio 6 and 20 cm but with a resolu- 
tion of 1" IHo fc UlvestadI i200lll detected a larger s cale com- 
ponen t extending to SW and NE from the nucleus. IVeilleu:j 
il99lll reported blue wings in the forbidden optical emis- 
sion lines and proposed a model with outflow and obscuring 



sured from the acquisition images (columns 5/6). 


Galaxy 


Morph. 


Seyfert 


Scale 


Image quality^ 




type 


type 


Pc/(") 


(") 


pc 


NGC 2273 


SBa 


2 


120 


1.00 


120 


NGC 3227 


SAB pec 


1.5 


84 


0.80 


67 


NGC 3516 


SO 


1.5 


183 


0.98 


180 


NGC 4051 


SABbc 


1 


45 


0.50 


23 


NGC 4593 


SBb 


1 


174 


0.49 


85 


NGC 4941 


Sab 


2 


72 


0.40 


29 



Although this has been obtained from the PSF measured in the 
acquisition image, we verified that the nuclear PSF measured in 
the reconstructed image of the Seyfert 1 galaxies is the same for 
NGC 4051 and only ~ 10% larger for NGC 4593. 

du st. The mid and far infrar e d flux es ha ve been explained 
by iRodriguez Espinosa et al.l lll996l) and IContini fc Viegai 
) ^ due to hot dust emission. In X-rays Lawren ce et al. 
1985r fo und variability on ti me scales of hundreds of sec- 
onds and ISalvati et alJ il993t) reported a flux change by a 
factor 2 at 2 .2 fim i n 6 m onths. Using X-ray maps of the nu- 
clear region ISinghl il999l) proposed the existence of a com- 
ponent associated with the nuclear activity and other asso- 
ciated with a starburst extended by ~ 40". 



2.5 NGC 4593 

This is a Seyfert 1 galaxy with detected nucle a r vari - 
ability in X-rays by IWeaver Gelbord fc Yaaoob J200l|) 
and in the optical and infrared b y [winkE^^L jj^^^, 
iKotilainen. Ward fc WilligeJ Jl993ft andlKotilainen fc Ware 
(1^^). The Ha map obtained by iGonzalez Delgado et al 
(|l993) shows ionized gas in a nuclear halo with radius 2-3" 
elongated N to S, and at the border of the halo some patches 
of emission which they suggest to be a broken starburst ring 
with minor axis length ~ 5" 7 running S to N. There are a 
number of H 11 regions in th e spiral arms but none in the 
large scale bar iEvans et alj[l996i'l . 



2.6 NGC 4941 



A ccording to the hard and soft X-rays spectra, 
iMaiol ino. Ris aliti fc Salv ati ( 1999) calculated that the nu- 
cleus of this galaxy has a column d ensity of Nh = 4.5 x 
10^^ cm~^. Extended [Olll] emission llPoggelll989l) shows an 
halo shape extending up to ~ 10" from the nucleus while 
Ha emission is concentrated in H 11 regions along the spiral 
arms. Schmitt ct al. (2001) found two resolved radio com- 
ponents separated by 15 pc along PA — 335°. 



3 OBSERVATIONS AND REDUCTIONS 

The data were obtained in queue mode over three 
semesters at the Gemini North telescope, using the GMOS 
IFU. The instrumental setup for each program is shown in 
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Table H The spectral resolution was R ~ 3000 (FWHM ~ 
100 km s~^) with a wavelength sampling of 0.692 A/pix. The 
GMOS IFU consists of an hexagonal array of 1000 lenslets 
which send the light through optical fibers to the spectro- 
graph. The lenslet array samples a field-of-view (FOV) of 
7" X 5". The sky is sampled with a field which is displaced 
by 1 arcmin from the object and has one half of the size of 
the object field. The centres of contiguous lenses are sepa- 
rated by 0'.'2, which is also the distance between opposite 
faces of the hexagonal lenses. The effective slit aperture for 
every lens is, however, O'.'Sl. 

Each half of the fibres (corresponding to half of the 
FOV) is aligned to a separate pseudo-slit. The data can be 
obtained using both slits (full FOV) or one slit (half the 
FOV, or 3'.' 5 X 5") with twice the wavelength coverage. The 
spectra are projected into an array of three 2048 x 4608 EEV 
chips disposed side by side with small gaps in between. The 
peaks of contiguous spectra are separated by ~ 5 pixels. 
In two-slit mode each pseudo-slit illuminates one side of the 
array, and the use of passband filters avoids spectral overlap. 

Before the spectral observations, one direct image was 
obtained for centering purposes, allowing an evaluation of 
the image quality. This value, measured as the FWHM of 
the stellar PSFs in the field, is listed in Tabled 

Data reduction was accomplished using generic IRAF^ 
tasks as well as specific tasks developed for GMOS data in 
the gemini . gmos package. The reduction process comprised 
trimming, bias subtraction, flat-fielding, cosmic rays clean- 
ing, alignment and interpolation of the data across the chips 
to recover all the spectra in one frame, extraction of the 
spectra (the tracing is done using flat spectra), wavelength 
calibration, sky subtraction and co-addition of different ex- 
posures. 

In Fig. Q we show, for each galaxy, the spectra from 
selected lenses. The left panels show the nuclear spectra, 
defined as the one with the strongest flux in the contin- 
uum around A8500. The right panel shows the spectra cor- 
responding to a lens located at l'.'7 from the centre of the 
array, towards the left side of the IFU field. 

The nuclear spectra of the Seyfert 1 galaxies NGC 4051 
and NGC 4593 present too much contamination from the 
AGN continuum and emission lines. In the case of NGC 4051 
these lines include the Ol blend at A8446, Ca-T in emis- 
sion and a number of Paschen emission lines. In the case 
of NGC 4593 the Ol emission is very broad and there are 
a number of other emission lines (e.g. Nl A8703.2). These 
strong emission lines precluded the measurements of the stel- 
lar kinematics within ~ 0'.'8 from the nucleus in these galax- 
ies. This effect was unexpect ed on the basis of th e previously 
inspected integrated spectra jNelson fc Whittl5 l995'l. prob- 
ably because it was diluted in their large aperture. In any 
case, the information we could obtain for these galaxies out- 
side the contaminated region is still valuable and have been 
included in the present work. The spectra of the Seyfert 1.5 
galaxies NGC 3227 and NGC 3516 (some authors classify the 
latter as Sy 1) present weak contamination from the broad 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agree- 
ment with the National Science Foundation. 



O I line which, however, does not affect the measurements 
significantly, as the two strongest absorption lines of the Ca- 
T are unaffected. 

Even after subtraction of the sky contribution, there re- 
mained some sky lines residuals in the extranuclear spectra, 
observed in Fig. Q For these lines as well as for the faint 
emission lines in the nuclear spectra of NGC 4941, we have 
set the parameters of the xcsao task (used in the velocity dis- 
persion measurements) to interpolate the continuum elim- 
inating all data points with ffuxes higher than 1 standard 
deviation of the average continuum flux within the cross- 
correlation window. We have checked that this procedure 
efficiently eliminated the residual emission lines. 



4 DATA ANALYSIS 

4.1 Velocity measurements; cross-correlation 
technique 

The kinematic measurements were performed by cross- 
correlating the spectra of individual lenslets with the spec- 
trum of a kinematic standard star observed with the same 
setup as the galaxies. The cross-correlation was performed 
using the task xcsao of the package rvsao iKurtz fc MinkI 
Il99di in IRAF, over the spectral range 8430-8900 A which 
includes the Ca-T features but avoids the noisy regions at 
the extremes of the spectra. The task uses the method of 
the quotient of the Fourier transforms applied to the spec- 
tra in velocity space and searches for the peak of the cross- 
correlation function, fitting a quartic function to the data 
with values above 50% of the peak and giving as outputs the 
peak velocity and the full width at half maximum (FWHM) 
of the cross-correlation function. 

As the surface brightness of the galaxies decreases to- 
wards the borders of the field, in order to improve the signal- 
to-noise ratio at these locations, we replaced each individual 
spectrum by the average of itself and the 6 nearest spectra 
(which are closer than O'.'S). This replacement was done for 
the spectra beyond a radius 1" from the centroid of the 
8500A continuum brightness distribution. 

In order to obtain the velocity dispersion values from 
the measured line widths we have convolved the spectrum of 
the kinematic standard star with Gaussian curves of known 
FWHM - where FWHM = 2.35 cr, and a is the velocity dis- 
persion - to create a set of synthetic spectra. We then per- 
formed the cross correlation between these synthetic spectra 
and that of the kinematic standard, measuring the FWHM 
of the cross correlation function as we did for the galaxies. 
By plotting the measured widths against the known Gaus- 
sian widths we obtain a very tight linear relation which is 
then used to obtain the real FWHM from the measured 
FWHM for each galaxy. Finally, as the radial velocities are 
measured relative to the kinematic standard star, the radial 
velocity obtained for each galaxy was corrected by the ob- 
servatory motion relative to the local standard of rest and 
by the standard star radial velocity as determined from the 
shifts between the measured wavelengths of the Ca-T lines 
in the standard spectrum and their rest wavelengths. 
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Figure 1. Typical spectra of our sample. In the left panels we show the nuclear spectra of the lens with the strongest continuum flux 
(adopted as corresponding to the nucleus) and in the right panels the spectra of a lens displaced I'.'T left from the centre of the array. 
Each left-right spectra pair corresponds to the same galaxy and is identified by a label in the left panel. The spectra have been brought 
to rest frame and the main spectral features are identified by vertical dashed lines. 
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Table 2. Details of the observations. 



PROGID Date Observed Field size Exp Time Grating/ Spec, coverage 

(MM/YYYY) galaxies (") (s) Filter (A) 



GN-2002B-Q-15 


12/2002 


NGC 3227 
NGC 3516 
BD +31 22 14 


3.5x5 


1440 
1440 


R400/RG610.G0307 
R400/RG610.G0307 


6400 - 


■ 10400 


GN-2003A-Q-20 


03/2003 


NGC 4051 
NGC 4941 
BD +21 24 25 


7x5 


1800 
1800 


R400/CaT.G0309 
R400/CaT.G0309 


8250 


- 9450 


GN-2004A-Q-1 


02/2004 


NGC 2273 
NGC 4593 
BD +31 22 14 


7x5 


3600 
3600 


R400/CaT.G0309 
R400/CaT.G0309 


7700 


- 9500 



4.2 Error calculation 

The IRAF cross-correlat ion task xcsao calcula tes the 
error following the method ofjTon^^^^avil llQTO") , which 
is similar to that used by iNelson fc Whittle! (1995). This 
method uses the signal-to-noise ratio R and the FWHM w 
of the cross-correlation peak, assuming that the error in a 
is the same as in Vr'. 

error = Aa = AK.^g^^. (1) 

For the spectra where the Ca-T lines are clearly de- 
tected and are not contaminated by emission lines, we ob- 
tain errors in the range 4-17 km s~^ for all galaxies. However, 
there is a loose correlation between the velocity dispersion 
measurements and the error calculated by this formula, due 
to a dependence in w that i s not fully canceled by (1 + R), 
an effect also pointed out bv lNelson fc Whittle! lll995h . This 
led us to adopt a conservative approach regarding the error 
reported by the xcsao task. After inspecting some represen- 
tative spectra and the calculated error values, we adopted 
15kms~^ as the errors in our individual measurements for 
Vr and a. We note that this is actually an upper limit, which 
may overestimates the actual errors for some of the measure- 
ments, in particular the ones from the more central spectra. 



find the largest difference. This is illustrated in Fig.|5| where 
we compare the a maps obtained with both techniques for 
NGC 2273. Note that the low dispersion ring (see discussion 
in Section fS.lH is clearly present in both maps, the only dif- 
ference being that the values obtained with the direct pixel 
fitting method are larger by ~ 20 — 30 km s~^. 

In order to investigate if this difference could be due to a 
template mismatch, we have performed velocity dispersion 
measurements using both tech niques and a se t of 2 8 stel- 
lar templates from the atlas of !Cenarro et all l|20^). The 
stars were chosen to span a wide range of spectral types. 
For the cross-correlation technique the results did not vary 
much for template spectral types G, K and M. Neverthe- 
less, as we moved to earlier spectral types the results began 
to vary but by simple visual inspection of these template 
spectra we verified that they are bad matches to the galaxy 
spectra. In the case of the direct pixel fitting method we 
made two tests. First we modeled each spectrum with a 
combination of all 28 templates, letting their proportions 
vary f reely, thus accounting for vari ations in stellar popula- 
tions iGarcfa-Rissmann et al . ''2005"). Then we repeated the 
fit using only our observed template star. The difference is 
negligible (at most a few km s~^), demonstrating that our 
kinematic measurements are not seriously affected by tem- 
plate mismatch. 



4.3 Alternative velocity measurements; pixel 
fitting technique 

In order to investigate the dependence of the veloc- 
ity fields and in particular of the velocity dispersions a 
on the measurement technique we ha ve also carried out 
a kinematic analysis in pix el space iBarth et al.! !2002!: 
!Cappellari fc Emselleinl!2004!) . 

Our implementatio n of the direct-fitting-method is es- 
sentially the one used in !Garcfa-Rissmann et aP l!2005!) . the 
only difference being that we rectified all spectra to account 
for the different slopes of stellar and galaxy spectra. 

We have compared our cross-correlation results with re- 
sults obtained with this method. We find that the radial ve- 
locity measurements are essentially unaltered. We also find 
that the patterns in the velocity dispersion maps remain 
the same. We do find a systematic difference in the veloc- 
ity dispersion values obtained from the two methods but 
this difference varies from galaxy to galaxy. For example, for 
NGC 4051 there is no difference at all while for NGC 2273 we 



4.4 Modeling of the radial velocity field 

The radial velocity maps obtained from the measure- 
ments described above (see Fig. I3I8II show that a clear rota- 
tion pattern dominates the stellar kinematics. In order to ob- 
tain an analytical description of the radial velocity field, we 
have adopted a very simple approximation, assuming plane 
keplerian orbits for the stars, and representing the gravita- 
tional potential of the bulge by a Plummer potential: 



which depends only on the scale length a and bulge mass 
M, and where r is the radial distance in the plane of the 
galaxy. 

We assume that the stars are in orbits close to a plane 
P{i, tpo) characterized by its inclination relative to the plane 
of the sky (i) and the position angle (PA) of the line of nodes 
-i/jQ. Although we do not expect that this simple approxi- 
mation provides a robust determination of the bulge mass. 
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Figure 2. Comparison of a maps obtained witli different methods for NGC 2273. The map on the left was obtained using the cross- 
correlation method, the one to the right was obtained using a pixel space fitting technique. In both cases we obtain the same main 
structure: a low a partial ring at 2" from the nucleus. 



inclination and bulge scale length, we can obtain reliable 
determination of the systemic velocity Vs, PA of the line of 
nodes ^po and the position of the kinematical centre. 

Using the above potential, we calculate the rotation 
curve in the plane P{i,tpo) with coordinates of the origin 
{Xq,Yo). The 2D rotation curve is then projected onto the 
plane of the sky where the radial variable is now R, the angu- 
lar variable is tp and the correspondent scale length is A. The 
relations between r and R, and between a and A are: r — aR 



sin^(^-i/»o) 



The 



and a = aA, where a = ^ coa'^{'ili — tpo) -i cas'^(i) 

systemic velocity Vs is added to the model as a zero point 
offset. 

In summary, there are 7 free parameters to be deter- 
mined by fitting the non-linear model 



Vr = Vs 



i?2 GM 



sin (j) cos [tj) — ipo) 



V I (^C0S2 (V' - Vo) + ) 



3/4 



(3) 



to the data Vr, what has been done using a Levenberg- 
Marquardt algorithm. 

Before performing the fits, we have inspected the map 
of uncertainties in radial velocity generated by the cross- 
correlation routine, the individual spectra and the radial 
velocity maps in order to obtain an uncertainty thresh- 
old to discriminate between reliable data (or measurement) 
and noise. This threshold varied somewhat from galaxy to 
galaxy, ranging from 9 to 17kms^^. Only reliable measure- 
ments were used in the data analysis and model fitting. The 
measurements considered unreliable have been masked in 
the velocity maps shown in Figs. I3I8I 

The parameters obtained from the fits are shown in Ta- 
bleO The kinematical centre (Xo, Vo) was used to calculate 
AXo = Xo - and AFo = ill - where and are 



the coordinates of the centroid of the continuum brightness 
distribution. 



5 RESULTS 

In this section we present and discuss the kinematic 
measurements obtained from the spectra as well as the re- 
sults of the modeling. 

2D maps of the measurements, together with the acqui- 
sition image and residuals of the model fits for each galaxy, 
are presented in Figs. I3I8I Each figure comprises 6 panels, 
as described below. 

(i) The top left panel shows a large scale /-band image 
with the IFU field drawn as a white box and the peak of 
the brightness distribution marked by the black cross. The 
small inset at the lower right corner details the morphology 
of the nuclear region covered by the IFU field. 

(ii) The top right panel shows the measured radial veloc- 
ities Vr (red crosses) along a virtual slit placed along the 
line of nodes and with width C/.'G, together with the corre- 
sponding model values (blue circles), as a function of the 
projected distance from the nucleus'^ R. The error bar at 
the lower right corner represents the maximum uncertainty 
in the velocity values. 

(iii) The remaining panels show 2D maps of measured 
properties. The middle left panel shows the flux integrated 
over the wavelength range used in the cross-correlation mea- 
surements. 

(iv) The middle right panel shows the radial velocity map, 



^ Note that, due to the fact that the slit width includes more 
than one lens there may be more than one velocity value for each 
given radius at the slit, even for the model. 
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Table 3. Parameters derived from our modeling. 



Galaxy 










M 


A 


i 




(km/s)^ 


(°) 


(pc) 


(pc) 


(xIO^Mq) 


(kpc) 


n 


NGC 2273 


1836 


233.4 


-7 ± 12 


-20 ± 12 


2.4 


0.17 


51 


NGC 3227 


1201 


149.7 


-2 ± 9 


-10 ± 9 


1.4 


0.11 


43 


NGC 3516 


2698 


48.4 


-18 ± 18 


9 ± 18 


15 


0.28 


22 


NGC 4051 


718 


106.5 


4 ± 5 


8 ± 5 


0.084 


0.031 


31 


NGC 4593 


2531 


88.5 


1 ± 17 


4 ± 17 


18 


0.41 


25 


NGC 4941 


1161 


-0.2 


1 ± 7 


-14 ± 7 


3.6 


0.24 


24 



The upper limit for the error in Vs is about 15kms ^. 
^ AXq and AVq are measured relative to the continuum centroid (see text). 



with isovelocity contours superimposed; white contours cor- 
respond to the measurements, while grey contours corre- 
spond to the best fit model. The kinematic centre is identi- 
fied by the black cross, and the line of nodes is indicated by 
the black dashed line. 

(v) The bottom left panel shows the velocity dispersion 
a map. 

(vi) The bottom right panel shows the residuals obtained 
from the subtraction of the kinematic model from the radial 
velocity map of panel (iii) . 

For display purposes, the kinematic maps of Figs. 
have been interpolated over the original hexagonal array to 
generate a rectangular grid uniformly sampled, with pixel 
size O'.'l X C/.'l. The scale (shown in panel (in)) and orienta- 
tion (shown in panel (v)) are the same for all the maps for 
a given galaxy. 

Table |1| lists the velocity dispersions measured from the 
integrated nuclear spectra inside circular^ apertures with 
diameters (/.'G and 1" for those galaxies with reliable nu- 
clear data, as di scussed in Section 3. Com paring our values 
with those from lNelson fc Whittle! (^9^, we conclude that 
they are compatible given the u ncertainties and o ur sm aller 
apertures. The value given by iTerlevich erall Jl99dl for 
NGC 3516 is significantly larger than ours , but these authors 
also p resent higher values than those of iNelson fc Whittle! 
il995h for other galaxies in common. For NGC 4051 and 
NGC 4593 the central velocity dispersions are unreliable due 
to the strong contribution from the AGN continuum. 

We now discuss the results obtained for each galeixy. 

5.1 NGC 2273 

The radial velocity map (central right panel of Fig. |^ 
shows a rotation pattern which deviates from the classical 
spider diagram, evidenced by the fact that the observed kine- 
matic minor and major axes are not perpendicular to each 
other, indicating deviations from axial symmetry. The resid- 
uals map show indeed deviations from the model above the 
noise level, which seem to be co-spatial with the star-forming 
ring (see discussion below). 

The rotation curve (upper right panel of Fig. 1^1 along 

^ Actually the combined spectra from all hexagonal elements 
with centres within the circular aperture radii. 



Table 4. Comparison between our a measurements and and 
those from the literature for the central spectra. Column 2 and 
3 presents our measurements (see text); column 4 and 5 present 
the values from the literature and corresponding apertures. The 
upper limit for the error in a is ±15kms^^. 



cr (0'.'6) a (1") cr (NW95) Aperture 
(km/s) (km/s) (km/s) (") 



NGC 2273 


104 




105 


136 


± 


22 


1.5 


X 


1.4 


NGC 3227 


114 




109 


128 


± 


13 


1.5 


X 


2.1 


NGC 3516 


192 




178 


235' 


i 




2.1 


NGC 4051 




85'' 




88 


± 


13 


1.5 


X 


2.1 


NGC 4593 




105'' 




124 


± 


28 


1.5 


X 


3.5 


NGC 4941 


138 




136 


109 


± 


13 


1.5 


X 


2.1 



From ITerlevich et al.l <199nl) . 

'' These values were obtained by inspection of the region around 
the contaminated central portion of the maps. See Sec. H^for 
details. 

the kinematic major axis seems to have reached the turnover 
at only 250 pc from the nucleus. In addition, the data 
turnover seems to be more pronounced than that of the 
model, suggesting a more concentrated mass distribution. 

The velocity dispersion map (bottom left panel of Fig. 
|3Jl presents a ring-like structure with values of about GO- 
SO km s~^ while in the centre the velocity dispersion reaches 
values of up to lOOkms"^. In Fig. |^ we present one- 
dimensional cuts of the velocity dispersion map along the 
major and minor axis of the galaxy which helps the visual- 
ization of the ring-like structure. This structure is co-spatial 
with the nuclear partial ring o bserved in the HST line emis- 
sion im ages and color maps of !Ferruit. Wilson fc Mulchaev! 
i2000l) . The lower velocity dispersion can be interpreted as 
a signature of recently formed stars, which still keep - at 
least partially - th e kinematics of the g as from which they 
have been formed JWozniak et alJ!2003!) and thus map the 
location of a starburst ring with semi-major axis of ~ 300 pc. 

5.2 NGC 3227 

Our radial velocity map (central right panel of Fig. |^ 
also shows a rotation pattern, where the redshifted side is 
to SE of th e nucleus, in agreement with th e H I ra dial veloc- 
ity data of !Schinnerer. Eckart fc Tacconi! l!2000!) . The line 
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Figure 3. Kinematic data for NGC 2273. The upper left panel shows a large scale i'-band image of the galaxy (spatial units are arcsec) 
with the IFU field overlaid as a white box and an inset detailing the nuclear region within the FOV of the IFU maps. The central 
left panel shows the flux integrated over the wavelength range of the cross-correlation measurements. The bottom left panel shows the 
velocity dispersion map. The upper right panel shows the measured radial velocities from a virtual slit placed along the line of nodes 
(red), together with the corresponding model values in (blue). The central right panel shows the radial velocity map with isovelocity 
contours from data and model superimposed (white for the measurements and grey for the model); the cross marks the kinematic centre 
and the dashed line shows the line of nodes. The bottom left panel shows the velocity dispersion map. The bottom right panel shows 
the residuals between the measured radial velocities and the model. All panels have the same orientation, indicated at the bottom right 
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Figure 4. Kinematic data for NGC 3227. Details are as described in the caption of Fig. |21 except for the fact that the IFU field is half 
the size, or 3'.'3 X 4'.'9. 



IFU stellar kinematics of active galaxies 11 




Figure 5. Kinematic data for NGC3516. Details are as described in the caption of Fig. |21 except for the fact that the IFU field is half 
the size, or 3'.'3 X 4'.'9. 
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Figure 6. Kinematic data for NGC4051. 



Details are as described in the caption of Fig. |21 
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Figure 8. Kinematic data for NGC4941. Details are as described in the caption of Fig. 1^ 
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Figure 9. One dimensional cuts of the 2D cr map of NGC 2273 
within pseudo slits CK'S wide crossing the nucleus at PA's 232° 
and 142°. 
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Figure 10. One dimensional cuts of the 2D a map of NGC 3227 
within pseudo slits 0'.'3 wide crossing the nucleus at PA = 156° 
and PA = 108° . 



of nodes PA obtained from the model (= 156.7°) agrees 
wit hin 4° with the values found by Mu ndell et al. ( 1993) 
and lSchinnerer. Eckart fc Taccoiiil l)200rfl . 

The residuals of the model fitting, shown in the bot- 
tom right panel of Fig. |11 are mostly below ±15kms^^, 
except for two blueshifted regions approximately along the 
line of nodes, which deviate from the model by more than 
20kms~^. The rotation curve (upper right panel) suggests 
that to NW the turnover has been reached already at 
~ 150 pc from the nucleus, while this does not happen to 
the SE. 

The velocity dispersion map (bottom left panel of 
Fig. 2J is characterized by a nuclear region with o above 
100 km s"'^ which drops in all directions but towards SE 
reaching values of « 75kms~^ or less. The a values increase 
again towards E, N and NW as evidenced in Fig. 1101 by two 
one-dimensional cuts of the cr map. The loci of low a re- 
gions have a good co rrespondence to the loci of high CO 
emission reported by ISchinnerer. Eckart fc Tacconil (|200(J) 
and referred to as a gas ring. In particular, the loci with the 
lowest velocity dispersions correspond to those with highest 
CO emission. The loci of low a i s also coincident with a ring 
of blue J — K colour reported bv lChapman et al] ll200(J) . Al- 
though our FOV seems to be too small to clearly reveal the 
rise of o beyond the ring in all directions, the coincidence of 
the low cr loci with both the CO and J — K rings supports 
also the existence of a low cr ring structure in our data. As 
in the case of NGC 2273 our interpretation is that we are 
observing a ring of recent star formation where the low cr of 
the stars is due to the low velocity dispersion of the cool gas 
from which they have formed. 



5.3 NGC 3516 

The radial velocity map (central right panel of Fig. I^J 
shows a rotation pattern which seems to reach the turnover 
at the edge of the field, corresponding to a distance of 500 pc 
from the nucleus, as can be also seen in the rotation curve 
(top right panel of Fig.jSJ. The velocity field is in good agree- 



ment with the 2D velocity map of lArribas et al.l il997ft . The 

kinematic major axis orientation of = 48.4° also agree 
with the values found by these authors. The above value for 
^0 is also the orientation of the photometric major axis ob- 
tained by fitting ellipses to the most external isophotes in 
the i'-band acquisition image (corresponding to an ellipse of 
semi- major axis 49"). 

From the rotation curve in the top right panel of Fig. |K| 
we obtain a peak-to-peak amplitude in the velocity field of 
^ 220kms~^. signifi cantly larger than the one measured by 
lArribas et alJ (1199711 (160kms~^). This discrepancy can be 
understood if we recall that those authors applied a l'.'5 
smoothing to the data, which would flatten the rotation 
curve. We note that the peak-to-peak amplitude of the stel- 
lar rotation is well below the one obtained from the gas 
kinematics (592 kms"^), bv lMulchaev et alJ (|l993). 

Our velocity dispersion map (bottom left panel of Fig. 
13 shows typical values in the central region larger than 
150kms~^ surrounded by a region with lower cr (as low 
as ~ 90km s~^). From measurem ents of the Mgib absorp- 
tion band, lArribas et al.l (^23) reported a central cr of 
164 ± 35kms~^ (within an aperture of 3"), which is con- 
sistent with our data, but a 2D distribution of line widths 
does not show any clear structure. In our data we can see 
lower cr values at locations surrounding the kinematic ma- 
jor axis direction at an average distance from the nucleus of 
~ 400 pc. 

The residuals from the radial velocity fit (bottom right 
panel of Fig. are mostly within the adopted error of 
±15kms~\ 



5.4 NGC 4051 

In this galaxy the nuclear continuum and broad Ca-T 
emission lines preclude reliable measurements of the stellar 
kinematics using the Ca-T absorption lines within 0'.'5 from 
the nucleus. This region has thus been masked for the model 
fitting. 

The radial velocity map (central right panel of Fig. |SJ 
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is dominated by the typical rotation pattern observed in all 
galaxies so far. The rotation curve (upper right panel of Fig. 

has a peak-to-peak amplitude of only 130kms~^ for an 
inclination of the galaxy i ~ 41° (assuming i — cos~^ b/a 
where a and b are the major and minor axis photometric 
lengths obtained from NED) and the turnover seems to oc- 
cur very close to the nucleus, at R ~50pc. 

The residuals from the model fit are within the 
rblOkms"^ error limits (bottom right panel of Fig. and 
do not show any well defined structure. 

The velocity dispersion map (bottom left panel of Fig. 
13, is also irregular with a values ranging from a minimum 
of ~ 60 kms~^ to a maximum of ~ 100 kms~^. Nevertheless 
it was not possible to extend the measurements closer to 
the nucleus than ~ O'.'S, making it difficult to draw any 
conclusion on the a distribution. 



5.5 NGC4593 

For this galaxy, as for NGC4051, the continuum from 
the Seyfert nucleus precludes reliable measurements of the 
stellar kinematics in the central region. 

The radial velocity map (central right panel of Fig. 
shows a clear rotation pattern in which the kinematical ma- 
jor axis runs approximately along E-W with the redshifted 
side to the E. 

The peak-to-peak amplitude of the velocity field is 
255 km s~^ , but the rotation curve (top right panel of Fig.[7|l 
suggests that the maximum velocity may occur beyond the 
edges of the observed FOV, which corresponds to ~ 700 pc 
from the nucleus along the kinematic major axis. 

The velocity dispersion map (bottom left panel of Fig. 
hints on a rise of the observed values towards the nucleus. 
At the most central but usable region we find values as high 
as 120kms~^. At ~ 2'.'5 from the nucleus (~ 400pc), the 
measured values drop to 80-90 km s~^, most notably along 
the kinematic major axis as in the case of NGC 3516. In the 
directions approximately perpendicular to the above the ve- 
locity dispersion is found to decrease slower, then rise again, 
reaching up to ISOkms"^ at the edges of the field. The map 
thus hints at the presence of a partial ring structure also 
for this galaxy. This is supported by the one-dimensional 
cuts along the major and minor axis of the velocity dis- 
persion map shown in Fig. 1111 Alternatively this region of 
low a values can be associated with a tightly wound spiral 
arm which can be observed in an HST F547M image from 
iMalkan. Goriian fc Tainl ^ll^^8^. 

iGonzalez Delgado et al have reported the pres- 

ence of a broken starburst ring in an Ha narrow band imag- 
ing. But a comparison of our data with theirs shows that 
our low a partial ring is nevertheless internal to their Ha 
ring. 

The radial velocity residuals (bottom right panel of Fig. 
|7|l are within the ± 15 km s~^ error limits, but show some re- 
semblance to the a map, suggesting systematic departures 
from circular rotation in association with the apparent par- 
tial ring of lower velocity dispersion. 

5.6 NGC 4941 



160 - 



140 



b 120 



100 





.PA = 


269° 


- W->E 


o 


oPA = 


20° - 


N->S 


o 








o 






o 


o 


• 

• 




o 

o o 



80 - 



-0.5 



CD a* • 
• o , 

On* 



9 • 



0.5 




kpc 

Figure 11. One dimensional cuts of the 2D cr map of NGC 4593 
within pseudo slits O'.'S wide crossing the nucleus at PA = 269° 
and PA = 20°. 



A rotation pattern is again observed in the stellar kine- 
matics (central right panel of Fig.jHJ. The rotation curve (top 
right panel of Fig.|HJ indicates that the maximum amplitude 
is located beyond the observed region, which corresponds to 
a radial distance from the nucleus of only 200 pc for this 
galaxy. 

The velocity dispersion map (bottom left panel of Fig. 
(HJ is almost fiat within ~ 100 pc from the nucleus, with 
central values ~ 135 kms~^, decreasing to ~ 100 kms~^ at 
the edges of the field (~ 200 pc). 

The radial velocity residuals (bottom right panel of Fig. 
ISJ are within the ± 15 kms~^ error limit and show no obvi- 
ous pattern. 



6 SUMMARY AND CONCLUSIONS 

In this work we have obtained 2D maps of the stellar 
kinematics of the inner few hundred parsecs of 6 nearby 
Seyfert galaxies at sub-arcsecond angular resolution, corre- 
sponding to spatial resolutions ranging from 30 to 180 par- 
secs at the galaxies. 

The stellar velocity field is dominated by rotation, well 
represented by a simple model where the stars follow plane 
circular orbits under a Plummer potential. The residuals 
between measured and modeled radial velocities are always 
smaller than ± 25 km s~^ . 

In most cases, the turnover of the rotation curve seems 
to occur within or close to the edges of the observed field, 
at radial distances ranging from ~ 50 pc for NGC 4051 
to ~ 500 - 700 pc for NGC 3516 and NGC 4593. Only for 
NGC 4941, for which our observations reach only as far as 
~ 200 pc from the nucleus, the turnover was not observed 
within the IFU FOV. 

The case of NGC 4051 is particularly interesting be- 
cause the turnover is at only ~ 50 pc from the nucleus, 
suggesting that the stellar motions are dominated by a 
highly concentrated gravitational potential. Indeed, the 
scale length of the Plummer potential obtained for this 
galaxy is the smallest of the the sample, A — 31 pc. 
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Adopting a = 85kms~^ for the nuclear velocity disper- 
sion (the mean value of innermost usable pixels) , and the 
iTremaine et alJ ||2002') relation, we infer a black hole mass 
of 4.5 X 10® Mq, which implies a sphere of influence of 
on ly 2.5 pc radius . Thes e values are consistent with those 
bv lPeterson et aD ||2q3)i who obtained a black hole mass 
of 1.9 X 10® Mq corresponding to a sphere of influence of 
1.2 pc radius. Such concentrated potential is unlikely to be 
responsible for the small potential-scale size we have found. 
A more likely interpretation is that the bulge is itself com- 
pact and concentrated. Indeed, unlike most Seyferts, this 
galaxy is late type (Scd), with a small bulge effective ra- 
dius. 

The velocity dispersion maps are featureless for 
NGC4051 and NGC4941. For the other galaxies the ve- 
locity dispersion maps show higher values at the centre and 
smaller values at 200 — 400 pc from the nucleus. In the case 
of NGC 3516 the lowest a values are observed towards the 
border of the IFU field, at location close to the kinematic 
major axis. For NGC 2273, the spatial distribution of lower 
0" values defines a ring-like morphology which is co-spatial 
with partial nuclear ring structures seen in line emission im- 
ages and color maps of previous works. For NGC 3227, the 
velocity dispersion map shows a very good correspondence 
between the region of low a values - which seems to delin- 
eates a partial ring structure - and regions of high ^'^CO 
(2-1) emission and blue J — K colours observed in previ- 
ous works. A partial ring of low a values is also hinted in 
NGC 4593 . As for the nuclear a-drops found by previous au- 
thors (e .g.[Emsellem et al.l200ll : lGarcfa-Lorenzo et aljl999l : 
iMarauez et alj|2003h the circumnuclear cr-drops we found in 
the present study can be interpreted as regions of higher 
(past or present) gas concentration, harboring younger stars 
which still preserve the lower velocity dispersion of the orig- 
inal gas from which they have formed. 
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